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We report on the development of a novel FPGA-based Time-to-Digital Converter and its imple¬ 
mentation in a detection chain that records the coordinates of single particles along three dimensions. 

The detector is composed of Micro-Channel Plates mounted on top of a cross delay line and con¬ 
nected to fast electronics. We demonstrate continuous recording of the timing signals from the cross 
delay line at rates up to 4.1 x 10® s“^ and three-dimensional reconstruction of the coordinates up to 
3.2 X 10® particles per second. From the imaging of a calibrated structure we measure the in-plane 
resolution of the detector to be 140(20) pm at a flux of 3 x 10® particles per second. In addition 
we analyze a method to estimate the resolution without placing any structure under vacuum, a 
signihcant practical improvement. While we use UV photons here, the results of this work apply to 
the detection of other kinds of particles. 


I. INTRODUCTION 

The detection of energetic particles, such as UV and 
X-ray photons or accelerated electrons and ions, is cen¬ 
tral to many research areas ranging from time-of-flight 
mass spectrometry [1] and diffraction experiments [2, 3] 
to the monitoring of metastable atoms [4] . A widespread 
technique to image these particles is provided by the use 
of a Micro-Channel Plate detector (MCP) which allows 
time-resolved two-dimensional imaging. In the above ap¬ 
plications both the resolution (in the two-dimensional im¬ 
age and along the time axis) and the maximum rate at 
which particles can be detected are of primary impor¬ 
tance. To obtain time-resolved two-dimensional images 
from the MCPs many methods have been developed [5]. 
In particular cross delay lines offer the possibility to reach 
high fluxes while maintaining a good spatial resolution 
along the spatial axes [6]. The capabilities of cross de¬ 
lay lines have been investigated in different contexts, in 
combination with MCPs [6-10], photocathodes [11] and 
Gas Electron Multipliers (GEM) [12, 13]. 

When recording three-dimensional information about 
the impact of incoming particles from delay-lines, the 
electronics which digitize and store timing signals are of 
fundamental importance. In particular, not only should 
the MCP and the delay line allow high fluxes but the 
entire detection chain, ending with the storing of parti¬ 
cle coordinates as digital data, must support these high 
rates. As the amount of data increases, both when lower¬ 
ing the time resolution to digitize the particle coordinates 
and when increasing the particle flux, the bottleneck in 
the detection chain is often related to the digital data. 
This limitation is either due to the time-to-digital con¬ 
verter (TDC) coding rate or due to the data transfer to 
a hard memory. In the works cited above, the maxi¬ 
mum rates at which data is recorded is a few 10^ s“^ 
[6-8, 13]. In [6] the authors anticipate read-out rates up 
to 1 X 10® s“^ but without a demonstration. In addition. 


while there exist several commercial TDCs digitizing at 
rates higher than 1 MHz per channel, we are not aware of 
any experimental evidence of a complete detection chain 
- ending with the reconstruction of 3D coordinates of in¬ 
coming particles - capable of acquiring data continuously 
at these rates. 

In this work, we report on the characterization of a 
complete detection chain to record the three-dimensional 
coordinates of incoming particles. The coordinates of 
a single particle are obtained by combining the timing 
information from the four outputs of the cross delay 
line. We have developed a FPGA-based TDC that al¬ 
lows us to demonstrate continuous acquisition at rates 
up to 4.1 X 10® s“^ per channel and to reconstruct 
three-dimensional coordinates of particles at rates up 
to 3.2 X 10® s“^, significantly higher than that previ¬ 
ously reported [6-8, 13]. An in-plane spatial resolution 
of 140(20) pm (averaged over the 80 mm diameter MCPs) 
has been measured at a flux of 3 x 10® particles per second 
by imaging a calibrated structure. This spatial resolution 
corresponds to a temporal resolution of 0.3 ns RMS on 
the measurement of time intervals. Building on the fact 
that the value of the resolution is not limited by the dis¬ 
cretization of time inherent to the use of the TDC, we 
discuss a quantity that allows us to estimate the resolu¬ 
tion in the absence of any calibrating structure. Using 
this quantity we investigate the performance of the de¬ 
tector at different incoming fluxes. We observe that the 
in-plane spatial resolution degrades with the rate of re¬ 
constructed particles [6]. 

This paper is organized as follows. We first present 
in section II the architecture and the performance of the 
novel FPGA-based TDC. The experimental set-up used 
to characterize the complete detection chain is described 
in section III. In section IV, we measure the count rates 
per output channel of the delay line and that of recon¬ 
structed particles with three-dimensional coordinates. In 
the last section V, we study the resolution of the detec¬ 
tor using two different approaches. First, a direct mea- 
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surement by imaging a calibrated structure is presented. 
Building upon this result, we introduce a quantity which 
yields, within the scope of our application, the spatial 
resolution and we investigate the influence of the flux of 
incoming particles on the measured resolution. 


II. OVERVIEW OF THE FPGA-BASED TDC 
A. Architecture of the FPGA-based TDC 

The time to digital converter we have developed is a 
FPGA-based TDC [14] using the Xilinx Virtex 4 technol¬ 
ogy. It features 17 independent input channels in which 
the arrival time of NIM signals can be digitized (see fig¬ 
ure 1). The digitization of time relies on a 16-stage delay 
line locked to a 520 MHz reference clock, so that the delay 
of each elementary step does not depend on temperature, 
supply nor aging (with no need of operating calibration). 
The duration of one elementary digitization step is equal 
to to = 120 ps, corresponding to one Least Significant Bit 
(LSB). Moreover, thanks to the locked delay line, the ref¬ 
erence frequency can be easily changed allowing one to 
modify the digitization step to (up to 1.9 ns). This per¬ 
mits one, for instance, to increase the range of the TDC 
or to distinguish distortions due to the TDC from those 
due to the detector itself. 

The time associated with the logical input signals 
(“stops”) entering the TDC is digitized over 26 bits in a 
numerical word of 32 bits. The duration of the recording 
can be extended from 8 ms, corresponding to the 26-bit 
period, to an arbitrarily long value with a mechanism 
permitting off-line counting of an elapsed time period. 
The words are stored in a FIFO register (named CHAN- 
nel REGister (CHAN. REG.)) with a depth of 509 words. 
Writing into the channel register results in an input dead¬ 
time of 2.5 ns. 

A coding sequence is enabled by an event trigger which 
is digitized in a similar word as the stops. The digitiza¬ 
tion of time is not deterministic because the TDC does 
not adjust the zero of its time base to the start of the ac¬ 
quisition. The TDC operates on a “sliding scale” mode: 
the time base runs freely and the time associated with 
a given stop is the result of the difference between the 
two time codes corresponding to the stop and start sig¬ 
nals. For an ideal TDC operating in sliding scale mode, 
the error in the determination of a duration is less or 
equal to 0.5 LSB, i.e. 60 ps RMS for the TDC we have 
developed. Conversely, the deviation of the mean value 
of many measurements of a given duration tends to zero 
with the number of measurements: in other words the 
TDC exhibits no “Integral Non Linearity (INL)” along 
its digitization range. Thus, the TDC does not induce 
any global distortion on the two-dimensional imaging in 
the MCP plane. 

The data words stored in the CHAN REG are multi¬ 
plexed at 16.4 MWord/s and stored alternatively in two 
acquisition FIFO registers with a depth of 509 words 



Output PCI bus 


FIG. 1. Architecture of the FPGA-based TDC. There are up 
to 17 input channels, each of which stores the time associated 
to an incoming NIM pulse into a CHANnel REGister of depth 
512 Words. The event trigger is time coded, stored and used 
as the reference from which time differences are evaluated off¬ 
line. The multiplexer is capable of treating 16.4 MWords/s 
and it stores the output data into two registers of depth 509 
Words (ACQuire REGisters). Data is then stored into the 
computer RAM via the PCI bus. 

(ACQ. REG.). The multiplexing is done by periodically 
reading the channel registers so that a data word is re¬ 
trieved from its channel register with the same probabil¬ 
ity independent of the channel. Finally, the acquisition 
registers are emptied at a 33 MWord/s rate and the data 
is transmitted through the PCI bus to the computer sys¬ 
tem. The stop times are digitally converted as long as 
the channel registers are not full: at high data rates, 
this implies that the number of converted stops from a 
given channel depends on the reading rate of the full data 
transfer. 

We have carefully optimized the depth of the ACQuire 
REGisters, on the one hand, and on the other hand, the 
PCI bus modes and wait states. This makes it possible 
to reach the maximum throughput of the multiplexer. 

The duration of the acquisition can be fixed externally 
by a start pulse (event trigger) and an end pulse sent 
to the TDC. In the data acquisition program, a polling 
loop runs continuously until all the acquisition registers 
of a run are transmitted to the computer RAM memory. 
The storage on the hard-drive disk and the reconstruction 
algorithm are executed when the acquisition is complete. 


B. Performance of the TDC 

To begin, we characterize the TDC by itself, indepen¬ 
dent of the rest of the detector. We first measure the 
Integral Non Linearity (INL ), the Differential Non Lin- 
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earity (DNL) and the resolution of the TDC. We check 
that the TDC exhibits no INL by sending pulses into a 
single input channel at a frequency of 20 MHz. The time 
interval coded by the TDC (50 ns) is stable to ±2 ps, 
a lower limit given by the instrument we use. We also 
investigate the stability of the INL as a function of the 
temperature. The temperature of the TDC is measured 
and stabilized within ±1 degree Celsius (10® samples). In 
a range from 35 to 86 degrees Celsius, the INL is stable 
within the precision ±2 ps of our measurement. 

In addition, we send pulse pairs with a fixed time in¬ 
terval into two separate channels to investigate relative 
drifts, both as a function of the temperature and as a 
function of the acquisition rate. At room temperature 
and for incoming rates varying from 0.4 to 4 MHz (see 
Fig. 2) we measure a drift of ^6 ps (on a time interval 
lasting 45 ns) negligible in comparison to the coding step 
to = 120 ps. 
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FIG. 2. Effect of the flux on the INL. The drift on the digi¬ 
tization of a fixed time interval (45 ns) between two separate 
channels of the TDC is plotted as a function of the repetition 
rate per channel. Inset: histogram of the counts on the TDC 
digit from a single measurement (flux of 4 MHz). 
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FIG. 3. Effect of temperature on the INL (a), DNL (b) and 
resolution (c) of the TDC. Temperature is measured and sta¬ 
bilized within ±1 degree Celsius (horizontal error bars). The 
vertical error bars reflect the level of accuracy of the measure¬ 
ment in the case of the INL (a) and the statistical deviation of 
the measurements in the case of the DNL (b) and resolution 
(c). 


In Fig. 3 a) we plot the evolution of the measured time 
interval (over 10® samples) and we observe a drift of ~ 10 
ps over the range 45 to 85 degree Celsius, corresponding 
to a small fraction (8%) of the coding step. 

The DNL quantifies the deviation of each digitization 
step from the ideal value; we measure a DNL (RMS 


value) less than 4% of the coding step to on each channel 
of the TDC (4 x 10® samples). In addition we monitor 
the variation of the DNL (10® samples) on a given chan¬ 
nel with temperature with the results plotted in Fig. 3 b. 
Over the temperature range we probe, the DNL is less 
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than 4%. 

The resolution cto of the TDC is defined as the stan¬ 
dard deviation of the distribution of repeated measures 
of a given duration. In the case of a non-deterministic 
TDC with a zero DNL uq is smaller than 0.5 LSB RMS 
[15]. It is well known that a non-zero DNL slightly deteri¬ 
orates the resolution [16]. We measure a mean resolution 
(To = 0.52(2) to over the 16 channels of the TDC, with a 
maximum value of 0.55 to RMS in a specific channel. In 
addition, the resolution (measured over 10® samples) is 
constant over a temperature range 40-85 Celsius degrees 
(see Fig. 3 c). 

To further characterize the TDC by itself, we measure 
the rate Pstops at which stops are recorded with an 80 
MHz digital pulse generator simulating the output pulses 
(NIM signals with 7 ns width) from the Constant Frac¬ 
tion Discriminator (CFD). The generator output is sent 
in parallel to 4 channels of the TDC to mimic the four 
outputs of the cross delay lines we use in the following, 
while the duration of the acquisition is determined by a 
second pulse generator. We observe 2 regimes: 

(i) when the amount of data does not exceed the 
depth of a CHAN. REG. (512 Words), we measure 
Pstops = 80 MWords/s per channel. An acquisition 
at this rate is limited to a duration of 6.4 ps. 

(ii) when the amount of data exceeds the depth of a 
channel register, we measure a rate Pstops = 4.12 
MWords/s per channel in agreement with the mul¬ 
tiplexer speed (16.4 MWords/s) over 4 channels. 

The role of the multiplexer is only visible in the latter 
regime. Whatever the regime, the dead time of a channel 
input is low enough not to play any role here. In the 
second regime the transfer rate through the PCI bus (33 
MWords/s) and the computing program polling speed are 
greater than the multiplexer speed. As a consequence the 
multiplexer sets the upper limit for the transfer rate of 
the TDC at Pstops = 16.4 MWords/s. 

III. DESCRIPTION OF THE TEST BENCH 

To characterize the detector as a whole we have 
set up a test bench consisting of an 80 mm diame¬ 
ter MCPs (Hamamatsu F1942-01 in chevron configura¬ 
tion) mounted above a cross delay line anode (Roentdek 
DLD80) and connected via vacuum feedthroughs to the 
electronics. 

We illuminate the MCPs (in vacuum at a pressure of 
2 X 10“^ mbar) with UV photons from a deuterium lamp. 
The emission spectrum of the UV lamp (185 nm to 450 
nm) is transmitted by a CaF 2 viewport (cut-off wave¬ 
length ^100 nm). The response of the MCPs strongly 
depends on the photon wavelength [17] with the conse¬ 
quence that the photons with the highest probability of 
detection have a wavelength A « 200 nm corresponding 
to 6-7 eV. UV photons are illuminating the entire surface 


of the MCP and the flux of photons is adjusted with a 
variable diameter hole situated roughly 850 mm from the 
MCP. The MCPs are biased with IkV voltage each and 
operate in the regime of gain saturation. 

Outside vacuum, the electronics are composed of the 
following elements. First, the four signals at the end 
of the delay lines are decoupled from the high voltage 
and amplified. The double stage amplifiers have 1 GHz 
bandwidth with a typical total gain of 20. Second, the 
signals are transformed into logical NIM signals through 
four Constant Fraction Discriminators (CFD) (Amplifier 
and CFD circuits are Surface Concept ACU 5-6-8.1). The 
CFDs have a dead-time of 7 ns. Finally the stop times are 
digitally converted with the FPGA-based TDC we have 
described above and are stored on the computer RAM. 


IV. CONTINUOUS ACQUISITION OF DATA AT 
HIGH FLUX 

We are first interested in characterizing the count rates 
at which the detector {i.e. the complete chain of detec¬ 
tion) works. We monitor two quantities, (i) the count 
rate per channel and (ii) the count rate of reconstructed 
particles. The count rate per channel is the number of 
stops per unit time recorded with the TDC on a given 
output channel of the cross delay line. The count rate 
of reconstructed particles is the number of particles per 
unit time with three-dimensional coordinates that one 
can correctly reconstruct based on the data acquired from 
the four delay-line channels. The event reconstruction is 
performed separately, and has no effect on the recon¬ 
structed count rate, therefore these computing resources 
are neither described nor characterized here. The mea¬ 
surement of the count rates is shown in Fig. 4 as a func¬ 
tion of the incoming flux of UV photons. We have plotted 
the maximum and minimum rates measured on the four 
channels with a blue area in between the extremes. The 
number of stops per channel is the same to within ±6% 
for each channel. The count rate per channel smoothly 
increases with the incoming flux of photons up to a max¬ 
imal value of ~ 4.1 X 10®/s in average. This maximum 
value is limited by the multiplexing task of the TDC, the 
speed of which is 16.4 MWords/s. 

The particle reconstruction is based on the selection 
of electronic pulses from the four channels and thus de¬ 
pends on the entire chain of detection as well as on the 
computer algorithm. The rate of reconstructed parti¬ 
cles we measure also smoothly increases with the photon 
flux up to a maximum value of 3.2 x 10® counts/s. As 
the photon flux increases, it deviates from the average 
count rate per channel. To quantify this behavior we 
plot in the inset of Fig. 4 the ratio rj of the rate of re¬ 
constructed particles to the rate of counts per channel. 
As it can be seen in Fig. 4, the count rate per channel 
slightly deviates from a linear increase at high flux. This 
is expected from the piling-up of pulses in the electronics 
and from the local gain depression of the MCP stack at 
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FIG. 4. Detected flux of counts (words) per output channel 
of the cross delay line (blue zone) and flux of reconstructed 
particles (red dots). The black dashed line corresponds to the 
maximum count rate per channel permitted by the TDC 16.4 
MWord/s transfer rate. Inset: efficiency rj of the reconstruc¬ 
tion process as a function of the incoming flux per channel 
(mean value of the four channels). 


high incoming flux [24] . The reconstruction rate strongly 
depends on these losses per channel as a reconstructed 
particle requires the digitization of all four pulses exit¬ 
ing the delay-line anode. In addition, the reconstruction 
algorithm we implemented considers all possible combi¬ 
nations of four digitized pulses and it does not induce 
any additional losses, a feature we checked. We thus at¬ 
tribute the decrease in the reconstruction rate at high 
flux to the observed increasing losses per channel. Nev¬ 
ertheless the fraction of the reconstructed particles lies 
above 80% over the entire dynamic range allowed by the 
TDC. We can compare these values with the performance 
in [7] where 50% losses were reported at a rate of 5 x 10® 
reconstructed particles. 

Importantly, the results in Fig. 4 were obtained while 
continuously acquiring data. With our long acquisition 
time (10 ms), a large number of transfers from the AC- 
Quire REGisters is necessary (typically > 150). This 
means that one can further extend the acquisition win¬ 
dow at no cost in count rate, allowing us to qualify this 
acquisition as “continuous”. The duration of the acqui¬ 
sition is not limited by the architecture of the TDC itself 
but by the size of the RAM as the digital words are con¬ 
tinuously stored in the computer. We have checked that 
measurements with longer acquisition times (100 ms, far 
above our final application needs), at the maximum rates 
yield identical results. In that later case, the number of 
transferred registers exceeds 3000. 

These high rates in continuous mode must be com¬ 
pared to the maximum achievable rate for the delay line 
with a 80 mm diameter MCP. It is well known that the 
detection efficiency of a MCP saturates when the output 
current becomes equal to a few percent of the strip cur¬ 
rent [18]. A rough estimate of the flux of incoming par¬ 


ticles at which the MCP we use saturates (5-10% of the 
strip current), yields typical values between 5 — 10 x 10® 
counts/s. Further improvement of the electronics may 
allow us to reach this limit. 


V. RESOLUTION OF THE DETECTOR 

In this section, we focus on determining the resolution 
of the entire chain of detection. We first present a direct 
measurement of the spatial in-plane resolution obtained 
by imaging a mask positioned slightly above the MCP. 
In the following section, we describe a measurement that 
allows one to evaluate the spatial resolution without the 
need for placing a mask above the detector. Throughout 
this section, we denote by x and y the spatial position of 
the incoming particle in the plane of the MCP surface. 
The position along the third axis is coded by the arrival 
time ta of the particle on the surface of the delay line. 

The recording of the timing pulses at the output of 
the cross delay lines, respectively tx, 2 ) for x and 
(ty,i, ty, 2 ) for y, allows one to reconstruct the three- 
dimensional coordinates of an incoming particle. The 
position X along the cc-axis is given by the time differ¬ 
ence (tx^l — ta:, 2)/2 with 

Ttot ^ 

ix,l — f-a 4“ ^ 

2 Vx 

rptot ^ 

tx,2 = ta - - -1-, (1) 

2 Vx 

where is the propagation delay from one end of the 
delay line to the other end and Vx = 1.05 mm/ns is an 
effective velocity describing propagation of the electronic 
pulses in the delay line [19]. Here the origin of spatial 
coordinates is set at the center of the delay lines. The 
position y along the j/-axis is extracted similarly from the 
time difference (typ—ty,2)/2. The coordinate ta along the 
third axis is determined from the time sums ^{txi +tx 2 ) 
or |(ty,i -I- tj/, 2 )- From these considerations, it is clear 
that the precision with which time is digitized is central 
to the resolution of the detector. 


A. Direct measurement of the in-plane spatial 
resolution 

A straight forward method to determine the spatial 
resolution in the x — y plane, i.e. the size of the point 
spread function in the plane of the detector, is to image 
an object with small structures onto the detector. We 
use a mesh mounted 4 mm above the front face of the 
first MCP plate. The mesh is a 0.05 mm thick Cu-Alloy 
sheet with a honeycomb pattern. The center to center 
separation is p = 0.95 mm with a wall thickness dp = 
0.08 mm. 

In figure 5 a), we show an image taken in the pres¬ 
ence of the mesh. In order to extract a spatial resolution 
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FIG. 5. a) Image of the honeycomb mesh; the blue line is 
the zone of data for the profile along x and the red line for 
the 60 degree profile analysis, b) Spatial prohle of the counts 
along the x axis, yielding a resolution Ux = 2.0(2) to. c) 
Spatial profile of the counts along the 60 degree axis, yielding 
a resolution creo = 2.5(2) to- 


from such an image, we plot the spatial profiles of the 
number of counts along two axes making an angle of 60 
degrees (see Fig. 5 b) and c)). We assume that the shape 
of the dark regions - corresponding to the walls - is lim¬ 
ited by the finite resolution of the detector. The profiles 
are fitted with a series of Gaussian functions, each with 
an identical RMS size and a separation given by the grid 
periodicity p [20]. We define the in-plane resolution of 
the detector along a given axis by this RMS width. We 
obtain ax = 2.0(2) to along the cc-axis and ago = 2.5(2) to 
along the axis making a 60 degree angle with the x-axis. 
The quantity to is the digitization step. In spatial units, 
this corresponds to 120(20) ^m and 150(20) ^m respec¬ 
tively. 


B. Estimating the resolution without imaging any 
structured object 

As discussed above, the time digitization of the sliding 
scale TDC sets a limit to the resolution with which a 
duration (evaluated as a time difference, such as At — 
tx,i — tx, 2 ) can be measured. The spatial resolutions Gx 
and Goo we have measured lie well above this lower limit 
Go ~ 0.5 to- Apart from the TDC itself, the sources 
contributing to the resolution are the electronic noise in 
the amplifiers and the walk (jitter) of the discrimination 
process. We assume all these sources of noise on a given 
channel can be considered independent from those of the 
other channels. We denote by at,noise the RMS width 
of the time distribution associated with a single channel 
(and therefore the measurement of a single time, such as 
tx,i). Since encoding a spatial position in the x — y plane 
involves time differences as defined in Eq. 1, we have 


The spatial resolution can thus be extracted from a 
measurement of Ct, noise- We introduce a quantity D to 
measure Ut,noise over the entire surface of the detector 
defined as 


D — (G,l + ix,2) — {ty,l + ty,2)- (3) 


D is derived from the four times measured from the de¬ 
lay lines. The RMS size gd of its distribution depends 
on the noise contributions from the four different detec¬ 
tion channels Ut,noise- Here, up to a constant geometrical 
factor, tx,i + tx ,2 — and ty^i + ty ^2 — are equal 
to the propagation time from one end of the delay line 
to the other end (see Eq. 1). For perfectly identical 
delay lines, cables and electronics, these two quantities 
are equal and D(x,y) = 0. In practice the situation is 
more complicated and the spatial distribution D(x,y) ex¬ 
hibits an offset and a gradient [21]. In order to extract 
a position-independent RMS size of the D distribution 
we proceed as follows. At a given position in the image 
plane {x, y), we record local histograms of the D{x, y) val¬ 
ues from which we extract a local mean value {D){x,y). 
The quantity tr^i is defined as the RMS size of the dis¬ 
tributions of D{x,y) — {D){x,y) summed over the entire 
surface of the detector. 

The direct measurement of the spatial resolution pre¬ 
sented in the previous section suggests that the contribu¬ 
tion of electronic noise is larger than that from the time 
digitization. Here we are interested in investigating these 
relative contributions to the quantity gd. To do so, we 
measure g^, as a function of the elementary digitization 
step Itdc of the TDC that we artificially change by re¬ 
moving one LSB at a time in the 32-bit word coding time. 
We have increased tTDC from its lower value tTDC = to 
up to ^tdc = 16 to- In Fig. 6 a) we plot the measured 
value Gd os a function of t^DC and for different fluxes of 
incoming particles. The pixel size, over which {D){x,y) 
is evaluated, is identical for all the txDC values (16 to)- 

We observe that gd increases as a function of txDC as 
expected since the digitization error increases with the 
coding step txDC of the TDC. In addition, a large offset 
value in gd is clearly visible at low txDC values, thus con¬ 
firming that the contribution of the uncorrelated noise 
<7t,noise dominates that from the digitization error. Fi¬ 
nally, we note that the quantity gd depends upon the 
incident flux: the larger the flux the larger gd (see be¬ 
low). 

The contributions to the RMS width gd from (i) the 
time digitization of the TDC and (ii) the uncorrelated 
sources of noise that we have defined previously, add in¬ 
dependently, leading us to write 

= O-g(txDc) + 4 CTj,noise (4) 


noise 


( 2 ) 


neglecting the contribution from the time digitization. 


where the first term corresponds to the error coming from 
the time digitization and the second term results from 
the independent sum of the contributions of the elec¬ 
tronic noise from the four input channels. Writing an 
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FIG. 6. a) Plot of an as a function of the TDC digitization 
step Itdc for three different fluxes of incoming particles. All 
quantities are normalized to the elementary step to = 120 
ps. b) Variation of an (circles) using the minimum digitizing 
step (Itdc = to) versus the incoming flux (expressed in re¬ 
constructed atom rate). The error bars reflect the variations 
of ctd as the pixel size is changed. The green diamond is the 
value of an obtained from the measured spatial resolution ago 
(Eq.5). 
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explicit expression for crg(tTDc) as a function of tg is not 
straightforward without further assumptions, but an ex¬ 
plicit expression is not required here. 

At low txDC values where the saturation of an is 
clearly visible (see Fig. 6a), the values of tTt.noise largely 
exceed that given by the digitization error, whatever the 
incoming flux of particles. This indicates that the resolu¬ 
tion is limited by independent noise in the different chan¬ 
nels of the detection chain. For the resolution Itdc = toj 
we conclude from our measurements (see Fig. 6) and 
Eq. 4 that ao — 2 Ut,noise- 

These considerations allow us to relate the quantity 
tTfl to the spatial resolution (see Eq. 2). We introduce a 


spatial resolution dres for the detector defined by 

_ ctd 

" 2V2' 


(5) 


The resolution tJres is a combination of the spatial resolu¬ 
tion along the two axis, and dy since, using Eq. 3, one 
can write cr|, = 4(tT^ -|- dy) and thus d^^^ = (tr^ -|- dy)/2. 
When the resolutions along the x- and y-axis are equal, 
dx — <7y, one obtains dres — <Jx — <yy 

At an incoming flux similar to that used in the 
direct measurement of the spatial resolution (200 — 
300 Kparticles/s), we find do ~ 7.7(4) tg and thus 
dres = 2.7(2) tg. To Compare this value to the direct 
measurement, obtained by imaging the grid, we evalu¬ 
ate dy in the following way. The resolution ugo along 
the axis making a 9 — 60 degrees angle with the x- 
axis reads dgg = (cos(9)dx)^ + (sin(9)dx)^- Using the 
measurement of dx and dgg we obtain dy = 2.7(1) to 

and \J{<7x + = 2.4(2) tg. This value is in good 

agreement with the spatial resolution Ures measured us¬ 
ing the quantity D, demonstrating that one can indeed 
use a measurement of do to evaluate the spatial resolu¬ 
tion from Eq. 5. The usefulness of the quantity do lies 
in the fact that it provides a good estimate of the spatial 
resolution without the need to install any object under 
vacuum above the MCPs or to otherwise illuminate the 
detector with a specific pattern. 

The contribution of the electronic noise (Tt^noise to the 
resolution along the third-axis coordinate ta can be eval¬ 
uated from the quantity D in the same manner as for the 
in-plane resolution. From Eq. 1 we obtain that the RMS 
size of the distribution of ta is dt^a = trres — 0.3 ns at 200 
Kparticles/s. We note that this low value might not set 
the resolution along the third axis ta in certain experi¬ 
mental situations with incoming particles impacting the 
MCP at low velocities [22] . 

Finally we investigate the influence of the incoming 
flux on the resolution do- The values of do are plotted 
in Fig. 6 b) as a function of the incoming flux of particles. 
One sees that do varies significantly over the range of 
particle flux studied: while it is constant in the lowest 
25%, its value is multiplied by a factor 2 at 80% of the 
maximum. A similar behavior has been observed with a 
lower threshold in a smaller (40mm) chevron piled pair of 
MCPs combined with wedge and strip anode [23]. This 
is related to the piling-up of pulses in the electronics and 
to the local gain depression of the MCP stack at high 
incoming flux [24]. 


VI. CONCLUSION 

In this paper, we have presented the characterization 
of a complete detector chain based on a pair of MCPs 
and cross delay-lines in conjunction with a FPGA-based 
TDC we have recently developed. This detector provides 
access to the three-dimensional coordinates of single par¬ 
ticles falling onto its surface at a high incoming rate. 










We have demonstrated that the architecture of our novel 
TDC permits a continuous acquisition of particles at 3.2 
million particles per second, a performance that is sig¬ 
nificantly better than previously reported to our knowl¬ 
edge (for similar space and time detector). A further 
improvement of this rate will be provided by future TDC 
electronic developments which aim at doubling the mul¬ 
tiplexer rate. 

We have measured the resolution of the detector by 
a traditional method (point spread function) and by us¬ 
ing a quantity D which yields similar results without the 
need to place any calibrated structure under vacuum. 
The simple use of D has allowed us to investigate the in¬ 
fluence of the detected flux of particles on the resolution. 


In addition, this simpler alternative can be important 
when it is impractical to introduce a grid as for instance 
when the detector is under ultrahigh vacuum [25]. 
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